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Abstract 

Old and new puzzles of cosmology are reexamined from the point of view of quantum the- 
ory of the universe developed here. It is shown that in proposed approach the difficulties 
of the standard cosmology do not arise. The theory predicts the observed dimensions of 
nonhomogeneities of matter density and the amplitude of fluctuations of the cosmic back- 
ground radiation temperature in the Universe and points to a new quantum mechanism 
of their origin. The large scale structure in the Universe is explained by the growth of 
nonhomogeneities which arise from primordial quantum fluctuations due to finite width of 
quasistationary states. The theory allows to obtain the value of the deceleration parameter 
which is in good agreement with the recent SNe la measurements. It explains the large 
value of entropy of the Universe and describes other parameters. 



1 Introduction 

The classical cosmology based on the equations of general relativity involving the principles 
of thermodynamics, hydrodynamics, the plasma theory and the field theory comes across 
a number of conceptual difficulties known as the problems of standard Big-bang cosmology 
[0, H, ^. These are the problems of singularity, size, age, flatness, total entropy and 
total mass of the Universe, large-scale structure, dark matter, isotropy of the cosmic 
microwave background radiation (cmb) and others. The various models were proposed 
for the solution of these problems. The inflationary model j^, ^ is the most popular 
one. There are alternative approaches which use the idea that in the early Universe the 
fundamental constants (velocity of light, gravitational constant, fine-structure constant) 
had the values different from the modern 10, Hi . 



1 



The observations of type la supernovae (SNe la) indicate that our Universe is accel- 
erating ^. This conclusion which appeared as partly unexpected for the cosmologists 
a few years ago nowadays practically does not called in question p||. The concept of a 
dark energy was proposed for the explanation of this phenomenon pv^ , 1 1 and the modern 
investigations in this field are directed toward filling of this idea with concrete contents 

The presence of the cosmological problems points to incompleteness of our knowledge 
about the Universe. It is generally accepted that the conclusions of classical theory of 
gravity cannot be extrapolated to the very early epoch. At the Planck scales one must 
take into account the quantum effects of both matter and gravitational fields. 

There cannot be any doubt that our Universe today contains the structural elements 
which bear the traces of comprehensive quantum processes in preceding epochs. The 
small cmb anisotropy and observed large-scale structure of the Universe Q can be given 
as necessary examples (see below). 

The application of basic ideas underlying quantum theory to a system of gravitational 
and matter fields runs into difficulties of a fundamental character which do not depend 
on the choice of a specific model. The problem of separation of true degrees of freedom 
under the construction of quantum gravity becomes of fundamental importance jist . It is 
commonly thought that the main reason behind such difficulties is that there is no natural 
way to define a spacetime event in general covariant theories [ [l^ . At the present time 
these difficulties are not overcome in the most advanced versions of quantum gravity. Also 
the quantum gravity cannot rely on experimental data . Therefore it is appropriate to 
construct the consistent quantum theory within the framework of the simple (toy) exactly 
soluble cosmological model. As it is well known the model of a homogeneous, isotropic 
universe (Friedmann-Robertson- Walker model) describes good enough the general prop- 
erties of our Universe. In this paper we study the model of quantum universe proposed in 
[ p^ p^ , pO| , PH . It docs not meet with the problems mentioned above and comes to the 
FRW model with positive spatial curvature within the limits of large quantum numbers. 

In Sect. 1^ we propose the method of removing ambiguities in specifying the time 
variable in the FRW model by means of modification of the action functional and find 
the solutions of obtained classical field equations. The Sect. ||is devoted to the quantum 
theory for a system of gravitational and matter fields. Here we formulate the equation 
which is an analog of the Schrodinger equation and turns into the Wheeler-DeWitt one 
for the minisuperspace model in special case. We concentrate our attention on study 
of the quantum universe which can be found in a region that is accessible to a classical 
motion inside the effective barrier formed by the interaction of the fields. We discuss the 
properties of the wave function of the universe and study the universe in low-lying and 
highly exited quasistationary states on the basis of exact solution of proposed quantum 
equation. In this section we calculate the proper dimension of the nonhomogeneities of 
matter density and the amplitude of fluctuations of the cmb temperature in highly exited 
state of the universe and propose a new possible quantum mechanism of their origin. The 
formation and evolution of large scale structure in the universe are considered as an effect 
of existence of primordial fluctuations due to flnite width of quasistationary states. The 
initially scale-invariant (flat) power spectrum of perturbations and the spectral index are 
calculated. The results are compared with the observed parameters of our Universe. The 
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flatness of the Universe and the large value of the entropy today receive their natural 
explanation. The observed accelerated expansion emerges as a macroscopic manifestation 
of the quantum nature of the Universe. 

Throughout the paper the notation Universe (with capital letter U) relates to our 
Universe, while universe (with small letter u) corresponds to arbitrary cosmological system 
of considered type. 



2 Classical description 

2.1 Coordinate condition and basic equations 

For simplicity wc restrict our study to the case of minimal coupling between geometry 
and the matter. Considering that scalar fields play a fundamental role both in quantum 
field theory and in the cosmology of the early Universe we assume that, originally, the 
Universe was filled with matter in the form of a scalar field cj) with some potential V{(j)). 
As we shall see the replacement of the entire set of actually existing massive fields by some 
averaged massive scalar field seems physically justified. We shall consider homogeneous 
and isotropic universe with positive spatial curvature. Assuming that the field (/) is uniform 
and the geometry is defined by the Robertson- Walker metric, we represent the action 
functional in the conventional form 



S = J dl] [iTa drjCJ. + 7r<j d,f(f) ~ H] . (1) 

Here r] is the time parameter that is related to the synchronous proper time t by the 
differential equation dt = N a dr/, where N{r]) is a function that specifies the time-reference 

scale, 0,(77) is a scale factor; tto and tt^ are the momenta canonically conjugate with the 
variables a and 0, respectively. The Hamiltonian H is given by 



E.= -N 
2 



■nl + ^nl-a'+a^Vic^) 



NU, (2) 

where the a is taken in units of the length I = ^/2/3^^lpl, I pi is the Planck length, and </> in 
units oi (j) = ■\/3/8 TT G. The energy density will be measured in imits = (9/16) rn'^pi- 

The function N plays the role of a Lagrange multiplier, and the variation 6S/SN 
leads to the constraint equation TZ = 0. The structure of the constraint is such that 
true dynamical degrees of freedom cannot be singled out explicitly. In the model being 
considered, this difiiculty is reflected in that the choice of the time variable is ambiguous 
(the problem of time). For the choice of the time coordinate to be unambiguous, the 
model must be supplemented with a coordinate condition. When the coordinate condition 
is added to the field equations, their solution can be found for chosen time variable. 
However, this method of removing ambiguities in specifying the time variable does not 
solve the problem of a quantum description. Therefore we shall use another approach 
and remove the above ambiguity with the aid of a coordinate condition imposed prior to 
varying the action functional. We will choose the coordinate condition in the form 

g^^(d^Tf = ^, or d^T = N, (3) 
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where T is the privileged time coordinate, and include it in the action functional with the 
aid of a Lagrange multiplier P 

S = Jdrj[TTad,,a + 7r^d,,<P + Pdr,T-n], (4) 

where 

n^N[p + n] (5) 

is the new Hamiltonian. The constraint equation reduces to the form 

p + n = o. (6) 

Parameter T can be used as an independent variable for the description of the evolution 
of the universe. Corresponding canonical equations reduce to the form 

2 

dra = -TTa, drTTa^ —Trl + a-2a^V{<j)), 



OtT = 1, drP^O. (7) 

Integrating the equation for P, we obtain P = E/2, where £' is a constant and the 
multiplier 1/2 is introduced for further convenience. The full set of equations for the 
model in question becomes |l^ 

{draf - y (drcl^f + U^E, (8) 

+ I (dra) (9t0) + ^ = 0, (9) 
a acp 

where U = ~ a'^'V {(/)). Equation represents the Einstein equation for the (g) 
component, while equation is the equation of motion Tq.^ — for cj), where TIf is the 
energy-momentum tensor of the scalar field 

n = Ti = Tl^-^{dT<Pf + V, 

Pi; = for /iT^z/. (10) 

From the analysis of the Einstein equations for this model it follows that inclusion of 
the coordinate condition (^) in the action functional leads to the origin of the additional 
energy-momentum tensor in these equations 

Tl^ = for fi^iy, (11) 
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that can be interpreted as the energy-momentum tensor of radiation. In the ordinary units 
E is measured in h. The choice of radiation as the matter reference frame is natural for 
the case in which relativistic matter (electromagnetic radiation, neutrino radiation, etc.) 
is dominant at the early stage of Universe evolution. If our Universe were described by the 
model specified by action functional (^) , it would be possible to relate the above radiation 
at the present era to the cmb. 

2.2 Solutions 

A feature peculiar to the model in question is that it involves a barrier in the variable a 
described by the function U. This barrier is formed by the interaction of the scalar and 
gravitational fields. It exists for any form of the positive definite scalar-field potential 
V{(p) and becomes impenetrable on the side of small a in the limit ^ — + 0. In general case 
{E ^ 0) there are two regions accessible to a classical motion: inside the barrier (a < ai) 
and outside the barrier (a > 02), where oi and 02 are the turning points (oi < 02) specified 
by the condition U = E. The set of equations (||) and (|^) determines the a and as the 
functions of time T at given V{(j)). When the rate at which scalar field changes is much 
smaller than the rate of universe evolution, i.e. {dtcj))"^ <C 2H^, where H = dta/a is the 
Hubble constant, and \df(j)\ ^ \dV/d(j)\, the equations (||) and (^) become 

{draf + U = e, (12) 



a ad) 



(13) 



where e and U depend parametrically on 0. In the zero-order approximation e = E. The 
solution to equation (H) can be refined by taking into account a slow variation of the field 
(j) with the aid of the equation 



y {dT<l>f + e{^) = E, 



(14) 



where e stands for a potential term. 

The solutions of the equation ([l3|) which determine the scalar field dynamics were 
studied in the inflationary models y, |^. The solution of the equation ([l^ ) at fixed value 
of (j) can be represented in the form 



ait) 



exp{-2Vy(t-i„)} 



1/2 



(15) 



where we denote 



(16) 
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Here a — a{tin) gives the initial condition for some instant of time t ~ tin- At a(0) = 
and a(tin) — 02 the corresponding scale factors are given in [|l8|, |l9j. The solution ( p^ 
shows that in the region a > a2 the universe expands in the de Sitter mode from the 
point a ~ a2, but in the region a < ai it evolves as a{t) ~ [2y/et]^^^ for 2y/V t ^ 1, 
that describes the evolution of the universe which density was dominated by radiation 
and as a{t) = ai — ^{t) with (^(t) ^ near the point of maximal expansion a = ai. The 
estimations for oi demonstrate that at small enough V the value a ^ ai can reach the 
modern values of the scale factor in our Universe. So, for the state of the universe with 
e ~ 1/AV and V ~ 10~^ GeV/cm'^ = 6.1 x 10~^^'^ (the mean matter-energy density in our 
Universe at the present era) we have ai ^ ^JXj^V ~ 10^^ ~ 10^* cm. 

In the extreme case of = 0, where there is no radiation, the region a < a\ contracts 
to the point a = 0, and the expansion can proceed only from the point a ~ a2 and the 
region a < a2 cannot be treated in terms of classical theory. Such models were widely 
enough studied by many authors (see e.g. ||, |2|, 23 



We concentrate our attention on study of the properties of the universe which is 
characterized by the nonzero values of E (and e) at the initial instant of time and can be 
found in a region that is accessible to a classical motion inside the barrier. 

The evolution of the universe depends on the initial distribution of the classical field 
and its subsequent behavior as a function of time. The solutions of the equation (^3[) 
for y ~ give evidence that the <\) decreases with time ^. From equations ^ and 
(p^, it follows that the inequality dxV + dT^/a^ < holds in the expanding universe. If 
V decreases with time, e can increase. Let us estimate e by using the relation e ~ T^a^. 
In our Universe, with a ~ 10^* cm, the main contribution to the radiation-energy density 
comes from the cmb with energy density p° ~ 10~^" GeV/cm^. Setting Tq = p2^, we find 
that, at the present era, the result is e = '--^ 10^^^ h. In the early Universe with a ~ 10~^^ 
cm and the Planck energy density we have e ~ fi-. It indicates that e should increase in 
the evolution process. This increase can be explained by a considerable redistribution of 
energy between the scalar field and radiation at the initial stage of Universe existence. 
Quantum theory is able to account for this phenomenon in a natural way as spontaneous 
transition from one quantum state of the universe to another. 

Taking into consideration the mechanism of quantum tunneling through the barrier 
and competing process of the reduction of V with time (which leads to the growth of the 
barrier U in width and height) allows to reexamine old and new puzzles of cosmology from 
the point of view of quantum theory. 

3 Quantum theory 

3.1 Quantization and properties of wave function 

In quantum theory, the constraint equation (^) comes to be a constraint on the wave 
function that describes the universe filled with a scalar field and radiation |jl^, |l^, ^ 

2zdT^=\dl~\dl~u\-^. (17) 
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Here the order parameter is assumed to be zero g2|, . This equation represents an 
analog of the Schrodinger equation with a Hamihonian independent of the time variable T. 
One can introduce a positive definite scalar product (^|^) < oo and specify the norm of a 
state. This makes it possible to define a Hilbert space of physical states and to construct 
quantum mechanics for model of the universe being considered. 

A solution to equation (Q7h can be represented in the integral form 



/oo 
dEei''^C{E)^E{a,(l)), (18) 
-oo 

where the function C{E) characterizes the E distribution of the states of the universe 
at the instant T = 0, while ijjE{o.,4') and E are, respectively, the eigenfunctions and the 
eigenvalues for the equation 

-dl + ^dl + U-E^^E = 0. (19) 

This equation turns into the famous Wheeler-DeWitt equation for the minisuperspace 
model H, |2|] in special case E — 0. 

A solution to equation (^9|) can be represented as 

/>oo 

V'B(a,(/.)= / de^,(a,<^) /,(</);£;), (20) 



where ipe and e are the eigenfunctions and the eigenvalues of the equation 

{~dl + U)^,^e^,. (21) 



For slow-roll potential V, when \d\nV/d(l)\^ <^ 1, the cp^ describes the universe in the 
adiabatic approximation and corresponds to continuum states at a fixed value of the field 
(j). The functions ipe can be normalized to the delta function S(e — e'). Their form greatly 
depends on the value of e. The quantities ff_{(j)]E) can be interpreted as the amplitudes 
of the probability that the universe is in the state with a given values of 4> and E . 

Since the potential U has the finite height Umax = and finite width then the 

quantum tunneling through the region oi < a < 02 of the potential barrier is possible. 
It results in that stationary states cannot be realized in the region a < ai. If, however, 
1^(0) <C 1, quasistationary states with lifetimes exceeding the Planck time can exist within 
the barrier. The positions e„ and widths r„ of such states are determined by the solutions 
to equation (^l|) for ip^ that satisfy the boundary condition in the form of a wave traveling 
toward greater values of a. Let us describe these states. 

We choose some value R > a2- Then ipf_[a,(j)) at fixed 4> can be represented in the 
form 

Lp^{a) = A{t) ^f\a) for < a < i?, (22) 

and 

ip,{a) = ^= \Lp[-\a) - S{e) ip^+\a)\ for a> R, (23) 
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where .4(e) and 5(e) are the amphtudes depending on e, (p^^^ is the solution that is regular 
at the point a = 0, normalized to unity, and weakly dependent on e, while ip[ \a) and 
ip)^\a) describe the wave "incident" upon the barrier (the contracting universe) and the 
"outgoing" wave (the expanding universe) respectively. Beyond the turning points the 
WKB approximation is valid so that one can write 

ip^^\a)^^ i cxpIt^ / VT^rfai^l. (24) 

The amplitude A{e) has a pole in the complex plane of e at e = e„ + jr„, and for a < i? 



the main contribution to the integral (20) over the interval — cxd < e < Umax comes from 
the values e ~ e„. The amplitude S{t) is an analog of the S-matrix 
The estimation 

/ 2 \'/' 

\^^-\a<R-[j^J W.X>R (25) 

shows that at r„ ^ 1 the wave function ^Pt{a) has a sharp peak for e = e„ and it is 
concentrated mainly in the region limited by the barrier. If e 7^ e„ then for the maximum 
value of the function we obtain 

where U (03) = and as 7^ 0. From this it follows that for r„ <C 1 the wave function reaches 
the great values on the boundary of the barrier, while under the barrier ~ 0(r„). 

In the limit of an impenetrable barrier, the function ipn = ipi'] reduces to the wave 
function of a stationary state with a definite value of e^. During the time interval 
AT < l/r„ the possibility that the state decays can be disregarded. This corresponds to 
defining a quasistationary state as that which takes the place of a stationary state when 
the probability of its decay becomes nonzero [2^. 



3.2 The universe in low-lying quasistationary states 

Calculation of the parameters e„, r„ of the quantum state of the universe can be done by 
both perturbation theory by considering the interaction a^V{4)) as a small perturbation 
against a? (in the region a < ai we have aP'V < 1) and direct integration of the equation 
( PH ) [ p8| , p^ . Such calculations show that the first level with eo = 2.62 — 1.31 fi. and 
To = 0.31 = 0.67 tp;^ emerges at V ^ 0.08 = 4.5 x 10"^ m%i. 

In the early universe, the quantity V{4i) specifies the vacuum energy density. The 
investigations within infiationary models suggest that the potential V{(j){t)) of the classical 
scalar field decreases with time. As the potential ¥{(/)) decreases, the number of quantum 
states in the prebarrier region increases but the decay probability decreases exponentially. 
The results of calculations are summarized in the table. 

Let us note that the quantum fluctuations of (j){t) in exponentially expanding universe 
can result in that the quantity (f){t) and the potential V ^ 0" will increase |^. Then 
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Table 1: Parameters e„ and r„ for various values of the potential V. At V = 
5.6 X 10"'^ rripi there are six levels in the system, three of them are displayed. 
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4.5 X 


10- 


-2 





1, 


.31 


6.7 X 10- 


-1 


2.8 X 


10- 


-2 





1, 


.40 


1.3 X 10- 


-2 


1.7 X 


10- 


-2 





1, 


.45 


4.3 X 10- 


-6 
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3, 


.17 


2.2 X 10- 


-2 


1.1 X 


10- 


-2 





1, 


.47 


1.5 X 10- 


-10 
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3, 


.30 


2.2 X 10- 


-6 
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4, 


.94 


6.5 X 10- 


-3 


5.6 X 


10" 


-3 





1, 


.49 


2.2 X 10- 


-24 
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3, 


.40 


2.2 X 10- 


-19 








2 


5, 


.26 


2.2 X 10- 


-15 



the quantum states of the universe in the prebarrier region cannot form. This case is not 
interesting for us and it wiU not be considered. 

The calculations demonstrate that the first instants of the existence of the universe 
(counted off from the moment of formation of first quasistationary state) are especially 
favorable for its tunneling through the potential barrier U . The emergence of new lev- 
els results in appearance of competition between the tunneling processes and transitions 
between the states. In the approximation of a slowly varying field 0, transitions in the 
system being studied can be considered as those that occur between the states \n) of an 
isotropic oscillator with zero orbital angular momentum which are induced by the inter- 
action a^V . In more strict approach which takes into account the variations of 1^(0) the 
transitions will be carried out at the expense of the gradient of the potential V{<j}) which 
follows from the quantization of the equation (|^) taking into account the evolution of the 
field 4> in approximation (^3|) [ p4[ . 

Considering the processes of transitions from some initial state m{Tf)) to final state 
n[T) (including the case m — n) and tunneling through the barrier from the final state 
as independent one can calculate the probability Wnm of transition between the states m 
and n as: 

W„™ « \{ipn\Ui\ip,n)\' exp{-r„ AT}, (27) 

where AT — T ~ Tq and Ui is the evolution operator in the interaction representation 
[ p6| . In the case of two- level system the computation of the total probability of universe 
decay, Wdec = 1 — (Woo + Ww), and the probability Ww demonstrates that over the time 
interval AT < 50 tp;, the transitions in the system predominate and only for AT ~ 10^ tpi 
the probability that the universe tunnels through the barrier becomes commensurate with 
the probability that it undergoes the — ^ I transition in the prebarrier region (see figure). 
Since the rate at which the level width r„ tends to zero is greater than the rate at 
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Figure 1 : Probabilities Wio and Wdec versus the time interval AT = T — Tq taken 
in units of the Planck time at F = 1.7 x 10~^ "^p;- 

which the V decreases, its reduction with time results in that transitions become much 
more probable than tunnel decays, in which case the former fully determine the evolution of 
the quantum universe in the prcbarricr region. If the universe has not tunneled through the 
barrier before the potential of the field decreases to a value V < 0.01 = 5.6 x 10~^ rup;, 
a sufficiently large number of levels such that the probabilities of decays from them can 
be neglected are formed in the system. Calculating the amplitudes of transitions over 
the time interval AT we find that the n n + 1 transition is more probable than the 
n — > n— 1, n+2 transitions. This means that the quantum universe can undergo transitions 
to ever higher levels with a nonzero probability. Since the expectation value of the scale 
factor a„ = (<Pn|a|(p„) ~ ^/n then it can be concluded that the characteristic size a„ of 
the universe that did not undergo a tunnel decay increases as it is excited to higher levels, 
i.e. the quantum imiversc being with respect to a in classically accessible region before 
the turning point ai can evolve so that a„ will increase with time. This can be interpreted 
as an expansion of the universe. The probability that, in course of time, the universe will 
occur in the region a > a2 outside the barrier is negligibly small. In the limit F — > 0, the 
universe is completely locked in the region within the barrier. 

The universe in the lowest state with n = has a "proper dimension" d sa ttoo ~ 
3 X lO"''"^ cm, the total matter-energy density p « 0.65 rripi, and the problem of the initial 
cosmological singularity does not emerge. The classical turning points are ai ~ 1.4 x 10"^'^ 
cm and 02 — 2.2 x 10~^^ cm. The value of ai determines the maximal dimension of the 
universe occurring in the lowest state in prebarrier region, while 0,2 characterizes its initial 
dimension after tunneling from this state. If a quantum universe tunnels from the states 
with n > 0, the dimensions of the region from which tunneling occurs can considerably 
exceed the Planck length. The constants E and V appearing in the Einstein equations 
will be determined by the corresponding quantum stage. 

Thus it turned out that quantum universe originally filled with a radiation and a 
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matter in form of a scalar field with a potential V{(j){t)) decreasing with time has nonzero 
probability to evolve remaining in prebarrier region. The expansion here is ensured by the 
transitions from lower states to higher states by means of interaction between gravitational 
and scalar fields. The system can be found in the highly excited states with n 1 as 
a result of such evolution. In states with n ^ 1 the potential V <^ 1 and r„ ~ 0. The 
state of the universe will be characterized by the quantum number n which determines its 
geometrical properties and new quantum number s responsible for the state of matter. 

3.3 Highly excited states 

The potential of will be chosen in the form V{(j)) — (m^/2) ip'^. From the condition y <C 1, 
it follows that the mass of the field must be constrained by the condition m <C mp; /|^|. 
In this case the states of the matter are determined by the solutions of the Schrodinger 
equation for a harmonic oscillator at given value of n . It describes the oscillations of (p 
near the minimum of the potential V{(j)). This process can be interpreted as the production 
of particles. At = 0, a similar mechanism leads to the production of particles by the 
inflaton field, which is identified with the scalar field Assuming as before that the 
V{4i) is slow-roll potential we find the condition of quantization of E 

E ^2N - {2Ny/^{2s + l)m, (28) 

where N — 2n + 1, and the values of the quantum number s are restricted by inequality 
s + 1/2 <C V2N^/m which reflects the fact that the mass m of the produced particles is 
finite. For small s the equality E « 27V holds to a high precision, so that the universe 
is dominated by radiation. A transition from the radiation-dominated universe to the 
universe where matter (in the form of particles produced by the field (j)) prevails occurs 
when the second term in ( ^8| ) becomes commensurate with the first one. The physical 
interpretation of the condition (|2|) will be considered below. 

For the universe with given quantum numbers n 3> 1 and s ^ 1 the wave function -ipE 
has the form |Q 

V's (»,</>) = ipn{a) f„s{(l)), (29) 

where 




This wave function is normalized to unity with allowance for the fact that the probability 
of finding the universe in the region a > 02 is negligibly small. 

The condition (EShcan be rewritten in the terms of "observable" quantities: the cosmic 
scale factor (a) = ^ N/2, where averaging was performed over the state (|2|) , and the total 
mass of the matter M = m{s + 1/2), 

E = 4(a) [(a) - M] . (30) 
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The classical universe is characterized by the total energy density p = Tq + Tq, where 
Tq and Tq are the energy-momentum tensors of the scalar field (^0|) and radiation (^Tj) 
respectively. Replacing all quantities by corresponding operators for the quantum universe 
we set ptot = (p)- It gives ptot = Psub + Prad, where 

193 M ^ E 

P^-' = l2W ^™'^^(^- ^^^^ 

Here in accordance with the Ehrenfest theorem we assume that expectation value (a) 
follows the laws of classical theory and the expectation values of the functions T" and Tq 
of a can be replaced by the functions of (a) . 

In the case, when p^^b 3> Prad we have (a) = M. This relation holds to a high 
precision ~ 10~^ in the observed part of our Universe, where (a) ^ 10^^ ^ 10^^ cm 
and M ~ 10^^ ~ 10^^ g. The quantum numbers of such universe are the following: 
n (a)^ ^ 10^^^ and s {a)/m ^ 10*° taking the proton mass for m. The value of n 
agrees with existing estimates for our Universe, while s is equal to equivalent number of 
baryons p7| . 

Thus for the matter-dominant era we have the following relation between (a) and Psub'- 



193 1 ^ ^'"^ 

12 Psub 



(32) 



On the other hand in accordance with assumptions made above in the universe with 
positive spatial curvature in this era the following equality must fulfill [^7| ^ 

(«> = (7y^— (33) 

\"0 ^ t Psub J 

where ilo — Psub/ Pc is the matter density in units of the critical density pc- From ( p2[ ) 
and (|3^ ) we find Q.q = 1.07. That is, the geometry of the universe with n ^ 1 and s ^ 1 
is close to Euclidean geometry (a flat universe). 

If one neglects the contribution from the kinetic term of the scalar field (tt^ = 0) 
then the corresponding ilo — 0.08. This value exceeds the contribution from the luminous 
matter (stars and associated material) [Q and it is close to the value for minimum amount 
of dark matter required to explain the flat rotation curves of spiral galaxies. Although the 
potential V{(j)) undergoes only small variations in response to changes in the field 0, the 
field (j) itself changes fast, oscillating about the point <^ = 0, so that the approximation 
in which tt^ = is invalid. The application of the present model in this approximation 
would result in the radiation-dominated universe; that is, it would not feature a mechanism 
capable of filling it with matter after a slow descent of the potential V{(i>) to the cquilibrimn 
position, which corresponds to the true vacuum. 

It is interesting to find the physical interpretation on (^). Passing on to the ordinary 
physical units we rewrite it in the form 

a^GMtou (34) 

where Mfot — M + Urad, U^ad = E/2a. It is easy to see that ( |3^ ) is the condition of the 
equality between proper gravitational energy of the thin spherical layer (with the total 



12 



mass Mtot) on the sphere with radius a and the sum of energies of particles M and energy 
of radiation Urad- In the modern era M ~ 10^*^ GeV 3> Urad ~ 10^^ GeV. If we extrapolate 
(^) on Planck era and set a — Ipi then Mtot = mpi. For the lowest quantum state we 
find that Urad ~ Tnpi/2. Since s = the parameter m « mpi. It means that vacuum 
energy of the scalar field and the energy of radiation make the comparable contribution 
to the total energy of universe with n = 0. In this state the scalar field \4>\ « 0.3 mpj. 



3.4 The nonhomogeneities of matter density 

The approach developed here makes it possible to obtain realistic estimates for the proper 
dimensions of nonhomogeneities of matter density, for the amplitude of fluctuations of the 
cmb temperature and points to a new possible mechanism of their origin, namely by means 
of finite values of the widths of quasistationary states. For a small, but finite value of the 
width r the quasistationary state does not possess a definite value of e. The corresponding 
uncertainty can serve as source of fluctuations of the metric 6a |20| . By associating 
e + (5e with the scale factor a + 5a and by using the solution for a(0) = we find the 
amplitude of fluctuations of the scale factor in the form 

5a 1 5ele ^^^^ 



a 4 1 -tanh(Vyt)/2VV^ 

Since 5e < F, the fluctuations 5a that were generated at the early stage of the evolution 
of the Universe will take the greatest values. For the lowest quasistationary state with 
e = 2.62, 5e < 0.31, V = 0.08, at i - 1 we obtain 

- < 0.04. (36) 

a 

Since the dimension of large-scale fluctuations changed in direct proportion a, this relation 
has remained valid up to the present time. For the current value of a ~ 10^^ cm we find 
that 5a < 130 Mpc. On the order of magnitude, the above value corresponds to the scale of 
superclusters of galaxies. Smaller values of 6e are peculiar to quantum states with smaller 
V. The fluctuations 5a corresponding to them are smaller than (^) and are expected 
to manifest themselves against the background of the large-scale structure. They can be 
associated with clusters of galaxies, galaxies themselves, and clusters of stars. 



3.5 Fluctuations of the cmb temperature 

The energy density of radiation can be expressed as 

Prarf = ^.g*T'', (37) 

where T is temperature and counts the total number of effectively massless degrees of 
freedom ^ Q . Using the relation (^) for prad we obtain 

E=^g.{aT)\ (38) 
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where we omit the brackets for simplicity. Leaving the main terms we can write 

ST 1 Se 6a , ^ 



For t <^ 1 it follows the estimation for the amplitude 

t Sa 

For the time t ~ 10^ yr corresponding to the instant of recombination of primary plasma 
(separation of radiation from matter), and for the observed value of e = 2.6 x 10^^^ h, for 
(H) we find 

— < 2.8 X 10"^ (41) 

Here VVt - 0.7 x 10"^. 

Upon recombination, the fluctuations of the temperature undergo no changes; there- 
fore, measurement of the quantity ST/T for the present era furnishes information about 
the Universe at the instant of last interaction of radiation with matter. The estimate 
in ( ^l|) is in good agreement with experimental data from which the trivial dipole term 



10 caused by the solar system motion was subtracted [29 



3.6 Large scale structure formation 

The problem of formation of structure in the Universe is nontrivial in any theoretical 
scheme In our approach the quantities (^6|) and (|4^) set a restriction from above 

on possible values of the amplitudes of fluctuations of cosmic scale factor and cmb tem- 
perature. In order to describe the power spectrum of density perturbations in the universe 
and the angular structure of the cmb anisotropics in the context of proposed approach it 
is necessary to have data about spatial distribution of the fluctuations 6e. For discussion 
we shall consider the mechanism of large scale structure formation in the universe based 
on the fluctuations Se which are distributed in space randomly. 

Let us assume that the perturbations Se depend on comoving space coordinates x = 



{x^,x'^,x^). Below we shall suppose that all perturbations "live" in flat space [31 and 
e-perturbations (5e(x) can be expanded in a Fourier series. According to Sect. 3.1 the 
state of the universe is characterized by the position e and width F of the level (here, for 
simplicity, the index n which specifies the number of level is omitted) . We assume that 
for given cosmic scale factor a the fluctuations (5e(x) have a form of Gaussian distribution 
in the coordinates {x^,x'^,x^) near fixed values Xq — {xq,Xq,Xq), 

where dispersions are supposed to be equal for three random values x^,x'^,x^. This 
distribution is normalized as follows 

'dx(5e(x)=r, (43) 
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where the integral is taken over space. Then the contrast (5c (x) = (5e(x)/e averaged over 
the whole space is 

r 

{S,{x))space ^ -■ (44) 

For averaged modulus-squared of the contrast 6^ (x) we have 

/li'k 
J2^nk), (45) 

where P{k) = |^e(fc)p is the power spectrum and ^e(fc) is the Fourier component of 5^{x.). 
For a homogeneous, isotropic universe 5e(fc) depends only on wavenumber k. From P{k) 
one can pass to the spectrum 



VM = ^P{k), (46) 



which is the measure of e-perturbations typical for the scale of wavelengths A — 27r/fc 
[|0[ ||. From (||) we obtain 

(|.-.(x)|"),,„.=ffnW=4,f^^. (47) 

Let us introduce the spectral index n{k) of the scalar e-perturbations as follows 

d\nP(k) 

The Taylor expansion of the spectral index about some fixed wavenumber ko 

n(fc)..(M+(^)jn- + ... (49) 
gives the following representation for the spectrum 



k 



n{ka) + {1/ 2) {dn/d In k)\ko ln{k/ko) + . 



P{k) = Piko) ^- j . (50) 

It shows that in a power-law approximation a scale-invariant Harrison-Zeldovich (HZ) 
spectrum Q corresponds to the case 

n{ko) = 1. (51) 

The spectrum P{k) can be expressed via the contrast (^^(x) 



P{k) = 



f sin(fca:) 



2 



(52) 
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Then the spectral index is 



i{k) = 2 



J dx cos(fca;) (5e(x) 



- 2. 



(53) 



Substituting ( [44 ) into ( ^2|) and (|5^) in the limit of small dispersions we find the following 
simple expressions for the spectrum and the spectral index 



P{k) 



sin(A:a;o) 



kxo 



n{k) = 2kxo \ cot(fca;o)| — 2, 



(54) 
(55) 



where = |xo|. As it is known [pT| , |3^ on the large scale the fundamental spectrum is 
consistent with HZ slope. In early epoch the equations ( pT| ) and (^5|) define primordial 
spectrum of e-perturbations with the wavelengths Xt = 2t: /k^ equal to 



Ai — 2.9x0, A2 = 1.4x0, A3 = 1.3x0, A4 = 0.82xoj 
A5 = 0.78 xo, Ae = 0.58x0, Ay^O.SGxo, etc. 



(56) 



In accordance with generally accepted views on mechanisms of formation of visible large 
scale structure in the Universe ^ one should choose the parameter xo equal to horizon 
which is determined by the width F as xo = l/F. Since the primordial spectrum is defined 
by the discrete set of wavenumbers k^ , then the HZ-spectrum itself can be written as sum 



over all possible roots (p6[) of transcendental equation (51), 



(57) 



Substituting (|5^) into (|£^) for the effective value of an amplitude of e-perturbations de- 
termined by root mean square of contrast 5^ (x) we obtain 



6"^ ^ 



1 



\/27r 



1/2 



(58) 



For the lowest state with e = 2.62, F — 0.31 and the wavelengths 
for the amplitude of e-perturbations of HZ-spectrum we find 



from and 



(59) 



The same estimation can be obtained if one makes a transition from integral in to 
sum over the vector k in a cubic lattice with spacing I/xq = F and then sums over the 
possible values of k for HZ-spectrum. 

The main contribution to ( ^9|) is made by the wavelengths A^ > xo. The amplitude 
(59) practically does not change up to the instant of recombination. Indeed, according to 
(54) we have 6^^ < (F/e) ^ 10""'^. Taking into account that amplitude of a-perturbations 
5a/ a remains constant during the evolution of the universe (see Sect. [3.4|) for the instant 
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of recombination from (B5|) we find that Sf^" < 10"^ and hence one can consider that the 
fluctuations ( p9[ ) also does not change up to the instant of separation of radiation from 
matter, so that 

Sf"'' - 10-2. (60) 

In order to relate the amplitude of e-perturbations with the density contrast 6p = dp/ p 
the expression for the energy density p = Tq + Tq we rewrite in the form 

P = (61) 

where we have used In the radiation-dominant era a ~ and 

Here we show in an explicit form the relation between our dimensionless and ordinary 
units. Since in this era to very high precision p ~ ||, H] the potential V in (^Tj) can be 
neglected. Then at given a for an effective value of energy density perturbations at the 
instant of recombination we obtain the following estimation 

After matter-radiation equality, the universe begins a matter-dominated phase and the 
density contrast dp increases according to known law 8p ^ (1 -|- z)-^, where z is redshift. 
As far as at the instant of recombination z = z^ec — 10'^ the perturbations (|63| ) guarantee 
obtaining the value (5p ^ 1 by now ^, |33|| . 

Thus in early Universe primordial e-perturbations which are distributed in space ran- 
domly can give the necessary value of energy density fluctuation during radiation domi- 
nation. Nonhomogeneities which arise here can grow up to observed large scale structures 
(galaxies and their clusters) in the Universe following the standard laws of general rela- 
tivity. At the same time there is no contradiction between the values (|6^) and (^l|). The 
amplitude of fluctuations (5T/T according to (|4^) takes into account the value of e which 
has the contribution only from "visible" cmb energy density. Whereas the value (^) ef- 
fectively includes the contribution from dark matter in the form of scalar field via the 
parameters e and F which are determined from the equation (pl|). 

Primordial fluctuations of (5e(x) present one of possible new mechanisms which can 
contribute to overall picture of formation of large scale structure in the Universe. 

It is interesting to clear up the possibility of the description of observed cmb anisotropy 
on the basis of the e-perturbations. This problem needs detailed study and we shall 
consider it elsewhere. In Appendix A we give some basic formulas in order to demonstrate 
in general the possible way of development of our ideas in this direction. 

3.7 Entropy 

The total entropy S per comoving volume [|l], ^ can be expressed as 



45 



g*. (aT)^ (64) 
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where g^,s can be replaced by for the most of the history of the Universe when all 
particles species had a common temperature. 

From ( |38| ) and ( |6^ it follows a simple relation between the E and the total entropy S 

For the adiabatic expansion, oT — const, the ratio E/S is conserved. Excluding aT from 
( |65| ) we obtain the relation. 

In the era with T ^ 10^^ GeV we have S ^ 1, but at present time for T ^ 10~^^ GeV the 
entropy S ^ 10^*. The large value of E today explains the large value of entropy of the 
Universe. 



3.8 Acceleration or deceleration? 

Recent measurements [0, ^ indicate that today the Universe is accelerating. Let us note 
that another possible explanation of the observed dimming of the type la supernovae at 
redshifts z ^ 0.5 is an unexpected supernova luminosity evolution At the present 

time the first interpretation of observed phenomenon is considered as more preferable. 

In terms of classical cosmology the accelerated expansion is described by the negative 
values of the deceleration parameter 

H'^ a 

In order to agree the experimental data with the theory it was proposed the concept of the 
dark energy which is nearly smoothly distributed in space. This dark energy component 
must have negative pressure that overcomes the gravitational self-attraction of matter and 
causes the accelerated expansion of the Universe. It is commonly assumed that the vacuum 
energy density in the form of a non-zero cosmological constant or due to a slow-roll scalar 
field called "quintessence" can be responsible for the dark energy jl^. [ij, [l3|] . 

Let us examine this problem from the point of view of the approach developed in this 
paper. To this end we rewrite (|67| ) in the form 

q = l + a — — . (68) 

K 

Bearing in mind canonical equation for SttTq from (|^) and having differentiated the equa- 
tion (|l^) with respect to a, we obtain that the derivative — drT^a must be substituted by 
the quantum mechanical operator 

n ^ i (-dl + ldl + a'~a^V- e) da. (69) 
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Then according to quantum mechanical principles the quantum analog of deceleration 
parameter can be calculated as 



('z) = i- 



(«n) 



(70) 



where averaging is performed over states "tpE and it is assumed that ofFdiagonal matrix 
elements from q vanish. (This corresponds to the representation of deceleration parameter 
as a scalar quantity.) 

In state with large quantum numbers n ^ 1 and s ^ 1, for the matter-dominant 
universe, where £'/4(a)^ <C 1, using the wave function (p9|) we obtain 



cos (27r(a)2) + -cos((27r-8) (o 
o 



(71) 



The expression in square brackets in (|7l| ) which contains two cosines rapidly oscillates 
with small period ^ "2,1 pi. Averaging ([7l|) over small interval near some fixed value of (a)^ 
we have 



(9) 



1. 



(72) 



This value can be associated with the deceleration parameter in classical theory. It agrees 
with the classical conceptions of general relativity about the expansion rate of the Universe 
in matter-dominated era with zero cosmological constant j2^, 28 1 . 

The quantity ( [72| ) does not take into account the quantum fluctuations of the scale 
factor 



(73) 



that specify root-mean-square deflection of the distribution |i/'£;(a, 0)p as function of a. 
In this case iI)e represents the wave packet which describes the universe being localized 
in space a near the expectation value (a) with deflection Aa. We shall show that at 
certain conditions (parameters of the universe) the fluctuations Aa can affect essentially 
the character of the expansion of the universe. It can provide in particular the accelerated 
expansion observed nowadays ||. 

We shall denote the scale factor taking into account the fluctuations Aa(t) as a(t), 
while the fluctuations themselves will be associated with quantity Aa(t) = a(i) — a{t), 
where a{t) is the scale factor without regard for fluctuations of considered type ([73|). 
Fixing some instant to for small intervals At = i — to we can write the expansion |27| 



a{t) = ao 



1 + ifo At - - go Hi + - So Hi At 



6 



(74) 



where sq = {1/ HQ){dfa/a)o and subscript indicates that corresponding values are taken 
at t = to- The similar series can be written for a(t) with the Hubble constant Ho, the 
deceleration parameter qq and so calculated with regard to fluctuations. It is natural to 
assume that the Hubble constant does not depend on the fluctuations (|73|), i.e. Hq — Hq. 
This assumption is based on astrophysical observations which do not record the necessity 
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to modify the classical conception of the Hubble's law. With regard to this facts from ( [74| 
and corresponding series for a(t) we obtain 

f Aa{t) Aao .Af^ , AaoV^ 1,^ Ntr2A^2 , 1 



ao ai J\ ao J 2 6 



Integrating (^5|) with respect to t from ti to to where \to — ti\ < Rc, Rc is radius of 
convergence of the series, we have 



\^ ao a§ " M " ' ao y 6 

- ^iso~so)Hi{to~hf + .... (76) 

Here a' and Aa* are the time average of the values a{t) and Aa{t) over the interval [ii, to]- 
Using the Einstein equations the parameter sq can be expressed in terms of qo, Qq and 
pressure po- Assuming that |f2o — i^ol ^ 1 and \pq — pq\ 'ti I, we find 

So - So = -{qo - qo). (77) 

For the instant of time when the universe is in the state with large quantum numbers, the 
mean (a) may be put to be equal to classical value ao- Then according to ( |73| ) it is natural 
to accept 

AH* = Aao = V(a2) - (a) 2 and a* w y. (78) 



Then up to discarded terms in ( fz^ ) 



go =90- , ,2 — fl + — ) fl-iHo(to-ti)) . (79) 

Using wave function (^9|) from ( |7^ ) we find that for the state with large quantum numbers 
the fluctuations 

^ = ^. (80) 

ao V3 

For such fluctuations 

For numerical estimation in the capacity of {to — ti) one can take the age of the Universe. 
For modern value to - ti = 14 Gyr [|| and iJo = GSkms'^ Mpc"^ we obtain 

qo = qo- 1.7. (82) 
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The parameter qo corresponds to the case when the fluctuations Aoq = and according 
to ( [r2| ) it equals to qo = (q) = I. As a result we find 

go - - 0.7. (83) 

This value takes into account the presence of quantum fluctuations of metric and it is in 
good agreement with SNe la observations |0, ^ . 

Let us note that for used values of Hubble constant and age of the Universe Ho {to — 
ti) < 1 and convergence of series ( [76| ) is not violated (see Appendix B). 

Thus, the observed accelerated expansion can be explained without implication of 
any additional concepts about matter-energy structure of the Universe considering this 
acceleration as macroscopic manifestation of its quantum nature. In any case unless the 
whole effect at least a part of it can be caused by quantum fluctuations of considered type. 

The represented calculations relate to the universe with large quantum numbers. In 
preceding epoch the Hubble constant and fluctuations Aa took different values. If one 
assumes, for example, that the relation H'^{t — ti)"^ ~ holds for earlier instants of time 
t, in epoch with h 1 the universe have to be decelerated if the fluctuations Aa < (a)/3. 
For more accurate calculation of the deceleration parameter of the universe in such states 
the averaging in ( |73| ) must be performed over the wave functions ipE which take into 
account that the variables a and in the equation (|l^) are not separated in general. 



4 Concluding remarks 

The main constructive element of our model which allows to avoid most of cosmological 
problems is an idea that E increased during the evolution of the Universe. The quantity E 
determines the energy-momentum tensor of radiation and can be found as an eigenvalue 
for the equation (p^. The above numerical estimations of the parameters of the quantum 
universe flUed with the radiation and scalar field show that the averaged massive scalar 
field used instead of the aggregate of real physical fields mainly correctly describes global 
characteristics of our Universe. It effectively includes visible baryon matter and dark 
matter. The kinetic energy term of the scalar field provides the modern value of the total 
energy density of the universe which is very close to the critical value. The status of 
the field (f> changes as we go over from one stage of universe evolution to another. In 
the early universe, the field (j) ensures a nonzero value of the vacuum-energy density due 
to V{(j)) values at which the equation ( |2l| ) for ip^{a,4>) admits nontrivial solutions in the 
form of quasistationary states. In a later era, when the field (j) descends to a minimum 
of the potential V{(j3) and begins to oscillate about this minimum, it appears to be a 
source of the particles of some averaged matter filling the visible volume of the universe, 
which has linear dimensions on the order of ~ (a). The galaxies, their clusters, and other 
structures in the Universe are subject to quantum fluctuations (due to the finite widths 
of the quasistationary states) that have grown considerably. 

The quantum fluctuations which specify the spread of the wave function of the universe 
in space of scale factor can ensure the accelerated expansion of the universe. In this 
sense they manifest themselves similar to dark energy. The theory gives the value of 
the deceleration parameter q = I (the universe is slowing down) for essentially classical 
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cosmological macrosystem and predicts q ^ —1 (the universe is speeding up) explaining 
the accelerated expansion as macroscopic manifestation of quantum nature of the universe. 
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Appendix A 

According to general approach (see e.g. ^) the detailed angular structure of the cmb 
anisotropy can be characterized by the two-point correlation function 

where is the angle between the directions ei and 62 in which the anisotropy is observed 
and the average goes over all points on the celestial sphere separated by an angle If 
one supposes that e-perturbations are distributed in space along the directions e, then 
according to (^5|), ( |39| ) for y^t ^ 1 the fluctuations of temperature in (Al) can be 
written in the form 

6T{e) I t 6e{e) 



T - 4 2V-e-t e ' ^^^^ 

From here it follows that correlation function C{'d) up to multiplier depending on time 
will be determined by 



;^^)o^=f:Q?PKcos^). (A3) 



1=2 

Here Pi is a Legendre polynomial, 

/ 

47r 



m— — l 

are the multipole moments and the coefficients aim are 



ai,n= ldn^-^Y*Je), (A5) 



where Y/m is a spherical harmonic, and the integral is taken over all directions in space. 
Specifying the form of distribution 6e{e)/e we can calculate the correlation function (j41). 
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Appendix B 



The first omitted term of the series (vWf has a form 



120 



(Fo-ro)ilo'(io-ii) , 



{Bl) 



where tq = {l/HQ){dfa/a)o and similarly for tq. In order to estimate it we shall suppose 
as in Sect. 3.8 that energy densities po, po and pressures po, po slightly differ from each 
other. Then we obtain that ratio 



TM (ro - ro) (to - 



^ (sq - So)H^ {to - ti) 
This value must be compared with ratio 

^ {so - so) Hiito ~ hf 



0.06. 



i {qo - qo) {to - h) 



0.23 



of first two terms of the series ([76|). These estimations show that the value qo 
be considered as reliable. 



{B2) 

(53) 
0.7 can 
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